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O ■ ABSTRACT 

We present the Spitzer IRS spectra for 33 young stars in Tr 37 and NGC 7160. The 
sample includes the high- and intermediate-mass stars with MIPS 24/im excess, the only 
known active accretor in the 12 Myr-old cluster NGC 7160, and 19 low-mass stars with 
disks in the 4 Myr-old cluster Tr 37. We examine the 10//m silicate feature, present in 
the whole sample of low-mass star and in 3 of the high- and intermediate-mass targets, 
and we find that PAH emission is detectable only in the Herbig Be star. We analyze 
the composition and size of the warm photospheric silicate grains by fitting the 10 fiva. 
silicate feature, and study the possible correlations between the silicate characteristics 



and the stellar and disk properties (age, SED slope, accretion rate, spectral type). We 



o 



find indications of dust settling with age and of the effect of turbulent enrichment of the 
disk atmosphere with large grains. Crystalline grains are only small contributors to the 



Q H i total silicate mass in all disks, and do not seem to correlate with any other property, 

except maybe binarity. We also observe that spectra with very weak silicate emission 
+3 , are at least 3 times more frequent among M stars than among earlier spectral types, 

which may be an evidence of inner disk evolution. Finally, we find that 5 of the high- 
and intermediate-mass stars have SEDs and IRS spectra consistent with debris disk 
models involving planet formation, which could indicate debris disk formation at ages 



as early as 4 Myr. 

' 



Subject headings: accretion disks — planetary systems: protoplanetary disks — stars: 
pre-main sequence 



1. Introduction 



The dissipation of protoplanetary disks around low- and intermediate-mass stars and the pro- 
cesses resulting in disk evolution and eventually planet formation are essential for the understanding 
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of the formation of our own Solar System. More than two decades of observations have allowed to 
pin down the main timescales for the presence of protoplanetary disks. Optically thick, accreting 
disks are nearly universal at ages around 1 Myr, become very rare at or after 10 Myr, and are 
absent in objects older than 30 Myr (Strom et al. 1989; Skrutskie et al. 1990; Lada et al. 2000; 
Carpenter 2000; Haisch et al. 2001; Armitage et al. 1999, 2003). The Spitzer Space Telescope 
has started to probe the processes of disk dissipation and evolution for large samples of targets 
(Porras et al. 2003; Meyer et al. 2004; Megeath et al. 2004) and to target the rare "transitional" 
disks, which are presumably evolving from optically thick disks to debris- like disk via dust settling 
and grain growth in the inner disk, and probably via the opening of physical inner gaps (see the 
examples of CoKu Tau4 in Taurus, Forrest et al. 2004; D'Alessio et al. 2005; TW Hya, Calvet et 
al. 2002; Uchida et al. 2004; among others). 

Observations suggest that disk evolution occurs faster in the innermost regions, as already 
predicted by Hayashi et al. (1985). The evidence of inner disk evolution is stronger in older 
regions, which also contain smaller disk fractions. In order to study large samples of evolved 
disks, we targeted the Cep OB2 Association (Patel et al. 1995, 1998). This bubble-shaped region, 
located at a distance of 900 pc (Contreras et al. 2002), contains the clusters Tr 37 and NGC 7160 
(Platais et al. 1998). Both clusters have been extensively studied with optical photometry and 
spectroscopy, and near- and mid-IR photometry (Spitzer IRAC and MIPS) in order to determine 
the cluster members and their characteristics (Sicilia-Aguilar et al. 2004, 2005,2006a,2006b, from 
now on Paper I, II, III and IV, respectively). 

Our previous work revealed more than 180 low-mass members in Tr 37 (spectral types G to 
M2) and over 60 in NGC 7160, with average ages 4 and 12 Myr and accretion disk fractions of 
~45% and ~4%, respectively (Paper II, III), belonging to the epoch when most disk dissipation 
and planet formation is believed to occur (Haisch et al. 2001). Combining observations at multiple 
wavelengths, interesting connections between the gas and dust evolution appear: grain growth 
and/or settling in the innermost disk affects nearly 95% of the objects, and the accretion rates fall 
in many cases well below the median rates in younger regions, to the point that some "transition 
objects" (objects with inner physical and/or opacity gaps resulting in no near-IR excesses at A < 
6/um, but optically thick disks at longer wavelengths; Paper III, IV) are consistent with no accretion 
(Paper IV). The accretion rates for the accreting transition objects are usually very low (<10 -9 
M /yr) compared to younger stars and to the bulk of accreting stars in Tr 37. This confirms their 
"transitional" stage, between Class II and Class III objects, and suggest a connection between gas 
and dust evolution. The transition objects represent about 10% of the total number of disks in 
Tr 37, and maybe a larger fraction at older stages, when disks are rare (Paper III, IV). The small 
fraction of "transition objects" around low-mass stars suggests a relatively rapid (up to ~ 10 6 yr) 
timescale for the dissipation of the outer disk after the inner disk is cleared. These short lifetimes 
may be related either to photoevaporation caused by the central star (Clarke et al. 2001; Alexander 
et al. 2006a,b; McCabe et al. 2006) or to (giant) planet formation at few AU (Lin &: Papaloizou 
1986, Bryden et al. 1999, Quillen et al. 2004). However, the presence of accretion and, therefore, 
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an inner gaseous disk in roughly half of these objects, suggests that photoevaporation alone cannot 
explain this behavior, at least not in most low-mass stars (Paper IV). 

Among the high- and intermediate-mass stars (spectral types BAF) identified in Tr 37 and NGC 
7160 (~60 in each cluster), our IRAC and MIPS data revealed the presence of a Herbig Be star in Tr 
37 and several candidates in both clusters, with 24 /im excess suggesting either debris disks (formed 
by second-generation dust after planetesimal or planet formation) or early-type "transitional disks" 
with large gaps. These disks would be among the few disks around BAF stars with well-determined 
ages (based on the low-mass cluster members) and, if they were confirmed to be debris disks, they 
would be among the youngest debris disks known. 

In order to study the characteristics of dust grains, the presence of grain growth and settling, 
and the nature of potential debris disks and "transition objects", we have targeted a sample of 
objects in Cep OB2 with the IRS spectrograph onboard Spitzer (see the list of targets and their 
characteristics in Table [1]). Our sample contains 20 low-mass stars (spectral types GKM) with 
different properties, and all the 13 high- and intermediate-mass stars (spectral types BAF) with 
any IRAC/MIPS mid-IR excess. In Section [2] we describe the data taking and reduction. In Section 
[3] we present the collection of IRS spectra and fit the silicate feature at 10/iin in order to estimate 
the grain composition and properties. In Section [4] we compare the IRS data with the available 
optical and IR information, and finally we summarize our conclusions. 

2. Observations and Data Reduction 

We observed the Cep OB2 objects as part of the GTO program "Spectroscopy of protostellar, 
protoplanetary and debris disks" in the IRS Campaigns 15, 16, and 22. They were observed with 
the nodded IRS Staring Mode. We utilized all three orders of the IRS Short-Low (SL) spectrograph 
(5.3-15 /im, A/AA ~ 90) and all three orders of the IRS Long-Low (LL) spectrograph (15-40 /im, 
A/AA ~90). This produces spectra covering the region from ~5.4 jum to ~35 /mi, which contains 
the most relevant dust features and traces the disk from ~0.1 to ~30 AU, approximately. In both 
clusters, we targeted all the high- and intermediate-mass stars with spectral types BAF and mid-IR 
excess (8 in Tr 37, 5 in NGC 7160), plus 19 low-mass stars in Tr 37 (about 1/8 of the total number 
of known members) and the only known accretor in the 12 Myr-old cluster NGC 7160. 

Given that parts of the two clusters do still show diffuse nebulosity that may result in some 
IR emission, we used the PRCS optical pickup on guide stars from the catalog for pointing. The 
exposure times were calculated individually based on the mid-IR fluxes measured with IRAC and 
MIPS. In order to be able to study the silicate feature at 10 /im, we required enough integration time 
to achieve a signal-to-noise ratio (S/N) of at least 30-50 in the 10 /im region. At longer wavelengths, 
our S/N requirements were only about 20, in order to optimize the GTO program. Except for a 
couple of stars, which had very low mid-IR fluxes, we achieved or surpassed our required S/N, being 
able not only to study the 10 /im feature, but to identify some other crystalline features in the 20 
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jum region as well (see Section 13. ip . 

We reduced the spectra using the Spectral Modeling, Analysis, and Reduction Tool (SMART; 
Higdon et al. 2004). We began with S12.0.2 BCD data from the SSC data reduction pipeline, which 
were already processed with stray light and flatfield corrections, and we departed from normal 
practice at several steps in order to achieve the best sensitivity and calibration. We corrected 
the pixels identified as NaN and rogue by interpolating the nearest non-bad neighbors (since, in 
case of bright and variable background, the correction of response in the rogue pixels by the SSC- 
recommended procedure of sky subtraction is not efficient). For each object that had multiple DCEs 
per exposure ID, the DCEs were averaged together. Next, we subtracted the sky by subtracting 
the off-nod position of the same order. We continued by extracting the spectra as a point-source 
for each nod position with a variable column extraction. We chose a width of 4.046 pixels at the 
center of each order and we let the widths vary from 3.2 to 4.9 pixels in SL2, 3.6 to 4.4 pixels in SL 
bonus order, 2.7 to 5.4 pixels in SL1, 3.2 to 4.9 pixels in LL2, 3.8 to 4.3 pixels in LL bonus order, 
and 2.5 to 5.5 pixels in LL1. This custom variable-width extraction window fits very tightly to the 
spectrograph point-spread function, so we can achieve good sensitivity. This point-source spectral 
extraction and calibration method was developed and refined for the IRS GTO surveys of Taurus, 
Chamaeleon and Ophiuchus. 

After extraction, the spectra of each nod were divided by the nodded spectra of our identically 
extracted photometric low resolution standard star a Lac, and then multiplied by the template 
spectrum of a Lac (Cohen et al. 2003). After reducing the data by this method, we merged the 
modules of each nod together, and then found the mean of the two nods. We determined the error 
in our spectra from the difference between nods. In a few of the objects, there are some small 
mismatches in flux between modules that are within our absolute spectrophotometric accuracy of 
10%. In these cases, we scaled the modules to match each other (see Table [2] for scaling factors and 
notes on the reduction). 

3. Analysis 

The IRS spectra of the Cep OB2 members are displayed in Figures [T] to [H For comparison 
and completeness, we display the available fluxes in the optical and near-IR (UBVRI, JHK from 
2MASS; Paper I, II), and IRAC/MIPS (from 3.6 to 24 ^m; Paper III). We also compare the 
resulting spectral energy distributions (SED) with the photospheric emission of a star with the 
same spectral type (Kenyon &: Hartmann 1995). It should be noted that all the IRS fluxes are in 
very good agreement with the IRAC/MIPS ones, finding significant deviations (over the 3-sigma 
level) only for the 24 fim flux of the star 01-580 in NGC 7160, for KUN-196 at 4.5 ^m, and for 
DG-39 at 4.5 and 8.0/mi. 

Figures [5] to [9] display a zoom in the 10/xm region of the spectrum. The silicate feature in 
emission is produced by warm (~300-500 K) dust grains of small sizes (~0.1-2 fim) in the optically 
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thin upper layer (disk atmosphere) of the disks (see Calvet et al. 1992; Chiang & Goldreich 1997; 
Menshchikov & Henning 1997, among others). The atmosphere of flared disks is heated by the 
stellar radiation, resulting in a higher temperature in the outer layers with respect to the midplane 
of the disk, which is responsible for the optically thin silicate emission in disks that are optically 
thick in the near- and mid-IR. Therefore, the absence of silicate features is associated to disks that 
have no warm small dust grains, either because the inner part of their disks (where irradiation can 
heat the disk atmosphere to the temperatures required) is evacuated from dust, and/or because 
there are no small dust grains in the disk atmosphere (due to grain growth and/or dust settling to 
the midplane). All the objects have spectra with good enough S/N at 10 /im to determine whether 
the silicate emission is present, although in a few cases the S/N and/or the strength of the silicate 
feature was not enough to allow proper fitting (see discussion in Section 13. 4p . The possibilities of 
finding silicate emission features at longer wavelengths (18 ^m dust resonance, or the crystalline 
features at 27.5 and 33.5 /im; see Jager et al. 1994; Henning et al. 1995; Waelkens et al. 2000; 
among others) is limited by our S/N in most of the stars. 



3.1. Disks, Dust, and SEDs in Low-Mass Stars 

Protoplanetary (or accretion) disks are optically thick, flared disks composed of gas and dust 
grains, generally extending from the innermost regions (~0.1 AU from the star) out to up to 
hundreds of AU. These are the disks typically observed around low- mass or classical T-Tauri stars 
(CTTS; Kenyon & Hartmann 1987; Chiang & Goldreich 1997, among others), and they can be 
characterized by excess at near- and mid-IR wavelengths and ongoing accretion. Nearly all the 
disks found around low-mass stars belong to this class, probably due to both relative numbers 
(timescale-dependent) and sensitivity limitations. As disks evolve, the IR emission from the inner 
disk tends to decrease (Paper III) and inner gaps may open and extend outward (Forrest et al. 
2004; Paper III), changing the disk structure so that the SED and the silicate emission are affected. 
Based on the presence of the silicate feature and the SED shape at different wavelengths, the disks 
around low-mass stars are usually classified as Class II protoplanetary disks (flared, with excesses 
at all wavelengths from near- to mid-IR, and silicate features in emission), transitional disks or 
"transition objects" (with no excesses at A <6^m, but CTTS-looking, optically thick disks at 
longer wavelengths, may have silicate emission), and debris disks (typically optically thin disks or 
rings, presumably formed by second-generation dust produced by collisions after planet formation 
has removed and/or conglomerated most of the disk material 0). 

A large excess over the photosphere at IRS wavelengths, together with excesses in the IRAS 
and 2MASS colors (see Paper III), characterize the CTTS-like disks. In Figures HE] and [5JH we see 
that the shape and strength of the silicate feature varies among stars with similar spectral types 



x Note that our limited sensitivity in our 24/im MIPS survey did not allow to detect such disks around stars later 
than F. 
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and similar SED slopes. This large variety of silicate emission had been noted in other studies (see 
Meeus et al. 2003; Kessler-Silacci et al. 2005). The fact that young, low-mass stars display such 
a large variety of disk morphology, dust size and dust composition, suggests that several factors, 
including probably age, spectral type, disk flaring, initial conditions, environment, and the presence 
of companions, affect the evolution of the disk. As examples, stars such as 11-2146, 11-2031, and the 
9 Myr-old 01-580 show very strong silicate features, characteristic of small amorphous grains; stars 
such as 11-2037, 23-162 and 14-183 show a weaker silicate emission, with some peaks associated 
to crystalline grains; and stars like 11-1209, 11-2322 and 11-2131 show nearly no silicate emission 
at 10 /mi over 3<r (see Section 13.41 for more details). In addition, some stars show amorphous 
and crystalline features at longer wavelengths, like 11-2131 and 11-2037 (both showing the 28/xm 
forsterite signature), 54-1547, 13-1250, and 21-2006. Note that the stars with detectable silicate 
features in the 18-30 /mi spectral range are among the higher-mass (G-early K) stars in our sample, 
which is probably due to the S/N limitations at A > 18 /mi for the faintest stars. 

Among the CTTS-like systems, our sample contains a double-lined spectroscopic binary (SB2), 
82-272, and two single-lined spectroscopic binaries (SB1), 14-183 and 12-2113 (Paper IV). The star 
82-272 shows remarkable emission features from crystalline and amorphous silicates, including the 
amorphous silicate feature at 18/mi and the crystalline forsterite bands at 23 and 28 /im. Meeus 
et al. (2003) suggested that the presence of companions may accelerate the evolution of the disk 
and the formation of crystalline silicates by inducing vertical mixing and stirring up the disk. 
Crystalline silicate emission at longer wavelengths is not detectable in 14-183 (maybe because of 
S/N limitations), but the 11.3/xm forsterite shoulder is marked (see Figured]). The third binary, 
12-2113, is similar to the rest of the stars in the sample. 

The IRS sample contains 3 of the transitional disks detected with IRAC/MIPS in Tr 37 (see 
Figure [T0j) . The star 14-11 is the most representative case of a non-accreting transition object in 
the IRS sample. Its IR excess starts in the region of the 10/mi silicate feature, very similar to CoKu 
Tau/4 (Forrest et al. 2004). Its age is about 1 Myr, again similar to the age of CoKu Tau/4 , as 
14-11 belongs to the young 1-Myr-old population associated to the Tr 37 globule (see Paper III for 
a discussion of the two populations in Tr 37). The silicate feature in 14-11 is very weak, being only 
visible due to the lack of continuum excess at shorter wavelengths. The IRS flux from ~10-30 /mi 
is roughly consistent with a black body with T~70 K, which would approximately correspond to a 
"wall" at a distance of ~10 AU for this Ml. 5 star. Given the lack of accretion signatures (so M< 
lO~ 12 M0/yr, Paper IV), it is a good candidate for hosting one or more large bodies or even planets 
(Quillen et al. 2004). 

The stars 73-758 and 24-515 are also transition objects, being the first a non-accretor, whereas 
the second is accreting at a very low rate (M ~ 10~ 9 -10~ 12 M /yr, marginally detectable in the Ha 
high-resolution spectra, Paper IV). The weak silicate feature in the transitional disks is consistent 
with a few-AU opacity gap, that would result in very little warm dust in small grains in case of low- 
mass stars. The K4.5 star 13-1250 had been classified originally as a potential accreting transition 
object (Paper III), although some small excess may be present at shorter wavelengths, suggesting 
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rather strong dust settling and/or grain growth as a mechanism to decrease the near-IR opacity 
rather than a physical gap. The presence of a relatively strong silicate feature in 13-1250 indicates 
that, despite the inner disk grain growth and settling that we infer from its small near-IR excess, 
there are still enough small dust grains in the disk atmosphere at distances ~l-3 AU. 

As mentioned before and explained in Paper IV, possible explanations of these gaps in tran- 
sitional objects could be either the strong dust coagulation up to planetesimals and maybe planet 
sizes (see Lin & Papaloizou 1986; Bryden et al. 2000; Quillen et al. 2004 among others), or the 
photoevaporation of the innermost disk by the central star (Clarke et al. 2001). Nevertheless, the 
presence of accretion (and hence, gas) in about 50% of the transitional disks (Paper IV) seems to 
favor strong coagulation to create the opacity gaps. The formation of one or more planet-sized 
bodies within the disk could lead to the opening of gaps that could be deep enough to prevent 
the flow of gas throughout the disk and would result in the end of the accretion processes as it 
presumably happens in CoKu Tau/4 (Quillen et al. 2004). Further IRS observations of the rest of 
transitional disks in Tr 37 (Paper IV) would be useful in order to determine the relation between 
the presence silicate feature (and thus, the presence of small grains in the inner disk) and other 
properties of the transition disks. 



3.2. Disks, Dust, and SEDs in High- and Intermediate-Mass Stars 

The high- and intermediate-mass stars (spectral types BAF) with disks are usually classified 
as Herbig Ae/Be stars (higher mass analogs of T Tauri disks; Henning et al. 1998; Dullemond et 
al. 2001; Bouwman et al. 2001; Meeus et al. 2003; van Boekel et al. 2003), transitional disks 
(with inner holes, also analogs to the low-mass transition objects), or debris disks (with inner 
gaps but basically depleted from gas and small dust particles after presumably planet formation, 
containing reprocessed second-generation dust replenished by collisions; e.g. Sylvester & Mannings 
2000; Rieke et al. 2005; Kenyon & Bromley 2005). In general, distinguishing optically thick, 
primordial, protoplanetary disks from optically thin, second-generation, debris disks at mid-IR 
wavelengths is not easy. Therefore, we will focus on the presence of silicate features as a proof of 
small grain content, characteristic of protoplanetary disks rather than of debris disks. Nevertheless, 
the weakness of the silicate feature in some of the optically thick disks around low-mass stars suggest 
that more detailed observations at longer wavelengths (including mm and submm), SED modeling 
and disk mass estimates will be required to probe the nature of the potential debris disks. Given 
the difficulties in distinguishing primordial and second-generation dust from mid-IR spectra only, 
we assign these disk the labels "transitional" or "debris" disks as the most probable nature of the 
objects, that will be studied in more detail with longer wavelength observations. 

As in the low-mass case, we observe very different spectral signatures and SED shapes in our 
sample of 13 high- and intermediate- mass stars (Figures [8] and [9]) . Here we list the characteristics 
of the most remarkable stars in Tr 37: 
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MVA-426: It is a B7 star in TV 37, found to have excess at all IR wavelengths (from JHK to 
70 fim) and strong emission lines in optical spectra (Contreras et al. 2002; Paper III). It shows the 
characteristic features of a Herbig Be star, with very strong continuum emission and remarkable 
PAH lines at 6.2, 7.7, 8.6, 11.3 and 12.7 /im, being the only star in the whole sample with detectable 
PAH features. It has no evidence of silicate emission at any wavelength, which could be either due 
to the reduced contrast with the high continuum level, or to the shielding of part of the disk by 
an inner wall. This makes MVA-426 very similar to HD 100453 (Meeus et al. 2003). Its IR excess 
emission is dominated by two "bumps" peaking at near-IR and mid-IR wavelengths, as it has been 
observed in other Herbig Ae/Be stars (see Hillenbrand et al. 1992; Meeus et al. 2003, among 
others). 

KUN-314s: It is the next star with the largest excess at 24 . This Ae star displays a 
featureless, relatively bright continuum spectrum, which is consistent with a classical Ae star with 
IR excess produced by the free-free emission of a gaseous disk (Porter & Rivinius 2003) or with 
a more evolved Herbig Ae star (Hillenbrand et al. 1992). It also shows an excess at all the IRS 
wavelengths and even a marginal 3-cr excess in JHK (Hernandez et al. 2006), which could suggest 
that it is a protoplanetary disk with strong sedimentation affecting the silicate emission. 

CCDM+5734Ae/w: This B3+B5 binary displays a weak but detectable silicate feature at 
10/xm, together with a relatively bright 18/xm silicate emission. The lack of an excess at wavelengths 
shorter than ~8 fim suggests a large inner gap in its non-accreting disk, and the presence of small 
silicate grains is consistent with a transition disk rather than with a debris disk. The fact that the 
18 iim. feature is comparatively stronger than the 10 emission, and that even the forsterite line 
at 33.5 ^m is marginally detectable, may suggest that the bulk of small silicate grains is located at 
larger distances from the star (in colder regions), probably due to the presence of the inner gap. 

KUN-196: This B9 star has excess at 5.8 and 8.0 pLva. IRAC channels, and a relatively smooth 
spectrum with very small silicate emission at 10 /im. This suggests the presence of a transitional 
protoplanetary disk with an inner gap and important dust settling and/or grain growth in its 
innermost region. 

MVA-1312: It is a B4 star with a featureless spectrum showing emission similar to (but lower 
than) the young debris disk of HR 4796 A (Jayawardhana et al. 1998), with no excess over the 
photosphere up to ~12 /im, and a smoothly increasing excess at longer wavelengths. The shape 
of the emission is compatible with a black body with T~100 K, that would mean a wall or ring 
at a distance of about 300 AU for a B4 star. Therefore, this object is our best candidate for a 
debris disk, although observations at longer wavelengths (mm and submm) are required to reveal 
the extent of the disk, the grain characteristics, and to find out if it is optically thin. 

The remaining intermediate-mass stars in Tr 37 (MVA-468, MVA-447 and BD+572356, spec- 
tral types B7, F0 and A4, respectively), do not show any silicate emission. The featureless spectra 
with small excesses starting at wavelengths from 8 to 20 /im suggest that they are most likely debris 
disks, analogous to the Vega-like systems described in Sylvester k, Mannings (2000). 
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In the 12 Myr old cluster NGC 7160, no high- and intermediate-mass stars were found to 
have optical emission lines (Paper II). The only star with a strong silicate emission, suggestive of a 
protoplanetary disk with a large inner gap, is the A7 star DG-481, which is the only one among the 
high- and intermediate-mass stars with a 10/im feature strong enough to be analyzed, (see Section 
ECU) . The stars DG-39, DG-682 and DG-912 (spectral types AO, A2 and F5.5) are consistent with 
objects surrounded by debris disks similar to the observed in Tr 37. The F5 star DG-62 does not 
seem to have a detectable excess larger than 3-cr over the photospheric level. 

3.3. The Shape and Strength of the Silicate Emission: Connections to the Inner 

Disk 

In order to analyze the silicate emission in a simple and systematic way, we normalize the 
silicate feature to the continuum level for all the high- and low-mass stars with significant excess 
at 10 /mi. The continuum is approximated by a straight line connecting the fluxes at 8 and 13 p 
(Bouwman et al. 2001), and the normalized silicate flux is defined as (Fy-F c )/F c , where "f" refers 
to the normalized flux in the feature, and "c" to the continuum normalized flux. This normalization 
helps to extract and compare the information contained in the silicate emission in a simple way, 
without assuming any particular model nor grain properties, and independently of the continuum 
level observed, and of the shape of the SED. Figure [TT1 displays the extracted normalized flux in the 
10 fim region for all the stars where the continuum can be successfully extracted (all the low-mass 
stars plus the A7 star DG-481). In order to examine the characteristics of the silicate emission, 
following the study of Herbig AeBe and T Tauri stars by Bouwman et al. (2001, 2006), we measure 
the peak of the normalized flux at different wavelengths (8.6, 9.8, and 11.3 /mi), corresponding to 
the broad shoulder noted by Bouwman et al. (2001), the peak of amorphous silicates with olivine 
stoichiometry, and the peak of crystalline forsterite. In addition, we measured the flux at the 
maximum of the feature. 

As a first step, we examine the flux ratios following the scheme designed by Bouwman et al. 
(2001, 2006). As stated there, we find a strong correlation between the flux ratios at 8.6, 9.8 and 
11.3 /im (see Figure [12]) . suggesting the same origin and similar composition for the silicate features. 
The values are different from the Herbig stars in Bouwman et al. (2001), but comparable to the 
ones observed for TTS in Bouwman et al. (2006) and in Kessler-Silacci et al. (2005). We do not see 
any difference in the behavior displayed by the A7 star DG-481, compared to the observations of 
T Tauri stars. We also observe a correlation between the silicate feature strength and the different 
flux ratios (Figure fT2"j) . similar to the correlation strength/shape described by van Boekel et al. 
(2003,2005), Przygodda et al. (2003), Bouwman et al. (2006), and Kessler-Silacci et al.(2005). 
The universal behavior of the relative fluxes, which do not depend on the spectral types nor other 
disk and stellar properties, seems to indicate that the components of the silicate feature and the 
processes leading to grain growth and dust evolution are very similar in stars with different spectral 
types (from A7 to M2.5). 
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Since the properties of the stars in our sample are well-known (see Table [T]) , we have explored 
the correlations between the different stellar and disk properties and the measured normalized 
fluxes. The significancy of these and other correlations investigated are listed in Table [5] There is 
a very strong correlation between the peak flux and the spectral type (see Figure [T2]) . which shows 
that the strength of the silicate feature decreases for the later spectral types. This effect could 
be related to differences in the region where the silicate feature is produced. Since late-type stars 
are less luminous, the areas where the irradiation is enough to heat the disk atmosphere to the 
temperatures required to produce the silicate emission are smaller and closer to the star for M-type 
stars compared to K-type and earlier. Therefore, the silicate emission in the lower-mass stars may 
be affected strongly by inner disk evolution (grain growth/dust settling/inner gap opening), or by 
the changes of the slope in the innermost disk (due to disk evolution) where the silicate emission is 
produced. Moreover, the size of the region producing the optically thin silicate feature is smaller 
for late-type stars, resulting in general in weaker emission. The A7 star DG-481 seems to be out 
of the general trend for the peak flux versus spectral type, but this is not surprising taking into 
account the difference in spectral type and in disk morphology with the rest of the sample. The 
silicate emission in an A7 star should be produced roughly at distances of 2-6 AU from the star, 
so relatively large gaps can develop before the silicate feature disappears. 

We also find a moderate correlation between the accretion rates and the silicate strength, 
showing that stars with strong features tend to have large accretion rates. The two stars with 
larger accretion rates in our sample (82-272 and 11-2146, with M~10~ M /yr) have very strong 
silicate emission. This is in agreement with the observations of FU Ori objects (Green et al. 2006), 
which are found to have strong silicate emission (although none of the objects examined here has 
accretion rates as high as FU Ori objects). The correlation between silicate strength and accretion 
rate may be related to the correlation between the spectral type and the accretion rate (Hartmann 
et al. 1998), given that stars with earlier spectral types tend to have stronger accretion rates (see 
Figure [13]) . We will explore this behavior in more detail in Section 14.11 

A possible dependence of the silicate feature on the SED shape, which affects the illumination 
of the disk by the central star, is explored in Figure [13] There is a general tendency for more flared 
(redder-colored) disks to have stronger silicate features, although we find strong scattering. This 
tendency can be understood as more flared disks are more efficient capturing the light from the 
star, so they can be heated to the temperatures required for producing silicate emission over larger 
areas. As noted in Paper III, we observe a general trend for a smaller flux at A <8 //m in Cep OB2 
compared to younger regions (i.e., Taurus) that could be interpreted as the effect of differential 
dust settling/grain growth (more efficient in the innermost disk) increasing with age. Given that it 
is in the innermost disk where most of the silicate emission for late-type stars is originated, if grain 
growth/settling/disk flattening is more efficient there, the silicate emission will tend to be weaker 
in later spectral types. 
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3.4. Fitting the 10 (im Silicate Feature: Grain Size and Mineralogy 



In order to analyze the 10//m silicate feature in a more quantitative way, we used the models 
developed for studying the objects in the FEPS sample (Bouwman et al. 2006), based on the 
previous work and observations of Bouwman et al. (2001) and van Boekel et al. (2003). These 
models, similar to the ones used by van Boekel et al. (2005), include 2 temperatures (for continuum 
and for grains), three different sizes for grains (0.1, 1.5 and 6.3 /im), and 6 different species (olivine, 
pyroxene, forsterite, enstatite, silica, and PAHs; see Table [3]). The models assume that the overall 
emission can be reproduced by adding the emission from the individual constituents so that the 
total flux can be written: 



We fit the spectral range from 6 to 13 /im, taking into account the errors at each wavelength 
and using a least-square minimization for obtaining the best fit. The fit could be extended to the 
silicate features in the 18-20 /jm region, but in our case the S/N in the LL part of the spectra 
is not high enough. There are some limitations to this approach. One of the limitations is the 
single temperature for the continuum and the single temperature for all the grains of different sizes 
and compositions, problem that could be improved by using detail disk models. The addition of 
the contributions of the different constituents is accurate if grains come in individual components 
or form very porous aggregates (see Bouwman et al. 2006 for details). In this study, we use 
this approximation as a first approach to characterize the grain sizes and composition, leaving the 
detailed disk modeling for future work. 

Spectra with very weak and/or noisy silicate features cannot be properly fitted. In our low- 
mass star sample, we find 4 stars with good S/N and non-detectable or marginally detectable 
silicate feature (11-1209, 11-2322, 11-2131, and 21-33), and 3 more stars with weak and noisy 
features (14-11, 24-515, and 21-2006), none of which could be properly fitted. Among the high- 
and intermediate-mass stars, we could only fit the silicate feature of DG-481, which is the only 
one showing a prominent feature at 10 /im. Figures [5] to [9] show the details of the silicate feature, 
including the fits whenever they were significant. Table S] summarizes the amount and size of the 
different materials needed to reproduce the silicate feature with this simplified models. It is worth 
mentioning that our fits reveal a very small amount of crystalline silicates (typically few percent, 
and never larger than 20% of the total mass, consistent with Meeus et al. 2003, Honda et al. 2003, 
2006, and Sargent et al. 2006), and no significant evidence for enstatite (only marginal detections 
within l-3a, consistent with the above mentioned references). 
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4. Discussion 

4.1. Parallel Dust and Gas Evolution: Grain growth, Turbulence, and Settling 

Using the grain properties obtained from model fitting (see Section 13.41 Table H|) , we examine 
the possible correlations between the characteristics of silicates, the age and the accretion rates. 
Figures [T4" l [151 [16l and [T7] show some of the correlations explored between the silicate size (obtained 
as a mass weighted average of the components required to reproduce the feature), the amount of 
mass in silicates, and the crystalline fraction. Their statistical significancy is compiled in Table [5j 
Even though the sample is small, since only 13 of the low-mass stars and 1 intermediate-mass star 
have enough S/N in the 10/mi region to allow size and composition analysis, we can draw some 
conclusions from the Tr 37 data. As already mentioned in Section 13.31 the behavior of the A7 star 
DG-481 is not significantly different from what we observe for the low-mass members, so we will 
analyze the whole sample together. 

There is a strong correlation that shows that the size of the grains tends to be smaller at older 
ages (Figure [14")) . This observation seems counterintuitive, since we would expect the grain size to 
increase with time as dust coagulates (Dullemond & Dominik 2005). Nevertheless, the 10 fim silicate 
feature traces only a (presumably small compared to the total) population of micron-sized grains 
in the upper layers of the disk. The population of grains that produces the silicate feature would 
be representative of the whole disk content only in the case of very strong turbulence (inconsistent 
with our moderate to low mass accretion rates) that would replenish the upper disk layers with 
coagulated and evolved grains formed in the deep layers of the disk. Therefore, this behavior is not 
inconsistent with grain growth to larger sizes occurring at later stages nor with inner disk removal 
(Paper III). Moreover, the lack of large grains (as traced by the silicate emission feature) at later 
evolutionary stages is consistent with dust settling increasing with time and affecting preferentially 
the larger grains. A KS test comparing the sizes of grains in disks with ages below 5 Myr and over 5 
Myr gives a probability of 0.07 (7%) that both distributions are identical. Taking into account the 
small size of the sample, the decrease of grain size with age is statistically significant. Since viscous 
disk evolutionary models (Hartmann et al. 1998) and observations of regions with different ages 
(Muzerolle et al. 2000, Paper II) show a decrease of the accretion rate with age, the decrease of 
silicate size with age would be a natural consequence if the evolution of the inner dusty disk occurs 
somehow parallel to the evolution of the gaseous component and turbulence/accretion (Paper IV). 

Moreover, the fact that we observe relatively large (nearly 6 (im) grains at very early stages 
may be an indication that grain growth (at least, to the sizes that we can still detect at 10 (im) 
starts already at very young ages (Dullemond & Dominik 2005). This observation implies that we 
must be careful before considering that stars showing a very "pristine" silicate feature (composed 
mostly of small grains) are very young or in an early evolutionary stage, since evolution (dust 
settling) may cause the large grains to be removed from the disk atmosphere, leaving only the 
smaller ones. A good example of this could be the only accreting star in NGC 7160, 01-580, which 
has a silicate feature due to very small amorphous grains, and an age near 9 Myr (5 Myr older 
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than the average star in TV 37). These observations seem in contradiction with the general idea 
that small grains are removed first in the atmospheres of Herbig AeBe stars because of coagulation 
(Kessler-Silacci et al. 2005), but ages of intermediate-mass stars should be regarded with care. 
Additionally, the average size of silicates is weakly connected to the slope of the SED (Figure [To]) . 
This result is similar to the observations in HAeBe stars in Bouwman et al. (2001), stating that 
the amount of silicate processing is not dominated by the SED shape, even if the overall strength 
of the feature may be affected by the disk flaring. 

Another strong correlation observed is that stars with higher accretion rates tend to have larger 
average silicate grains (Figure [T4"|) in the disk atmosphere, compared to stars with low accretion 
rates or even no accretion. Given the relationship between accretion and turbulence (Hartmann 
et al. 1998; Schrapler & Henning 2004), and that the silicate feature traces the grains in the 
warm disk atmosphere, these observations could be an indication that more turbulent disks can 
replenish the disk atmospheres with more (and larger) grains that might otherwise settle in the 
middle plane in absence of strong turbulence. Nevertheless, the scattering is large, and, as we 
mentioned before, some of the more massive stars (82-272, 11-2037) may have large accretion rates 
because of their more massive disks, and moreover, strong silicates due to heating of larger disk 
areas to the temperatures required for producing silicate emission. The fact that we observe more 
mass in warm silicates in stars with higher accretion rates (Figure [T6j) is likely to be due mostly 
to the higher accretion rates around earlier spectral type stars, since the total mass in silicates is 
mostly correlated to the size of the area producing the silicate feature. As we will discuss later 
(Section I4.2p . the stars with very weak silicate emission are preferentially M-type, so effects such 
as illumination or inner disk evolution may be more important for them. But, in any case, the 
correlation also holds for K and earlier spectral types, so for a given spectral type, the silicate 
emission may reflect the influence of turbulence in the grain content of the disk atmosphere. 

We do not observe any correlation between the crystalline grains or the degree of crystallinity, 
and any other property (age, accretion rate, SED slope, spectral type), except for the fact that the 
two objects with more than 15% of the mass in crystalline grains (Figure fl~7l) are the spectroscopic 
binaries 82-272 (SB2) and 14-183 (SB1). The double-lined spectroscopic binary 82-272 (Paper IV) 
has two spatially unresolved accreting components, and a very strong silicate feature at 10 (im, 
plus relatively strong amorphous and crystalline features at 18-28 fim. The single-lined binary 14- 
183 has strong crystalline features at 10 //m as well, but the third SB1, 12-2113, has a crystalline 
fraction similar to the rest of the stars, which have all less than 10% of mass in crystalline grains. 
As mentioned by Meeus et al. (2003), the higher crystallinity fraction in close binaries could be a 
consequence of their stronger activity and more turbulent disks. Nevertheless, the other effect of 
enhanced dust processing mentioned by Meeus et al. (2003), that binaries would have larger grains, 
is not observed here, although we must keep in mind that the silicate feature traces only the upper 
layers of the disk, and large grains can be removed from the disk atmosphere via dust settling. 

The lack of correlation between crystalline silicates and other disk and stellar properties has 
been noted by Meeus et al. (2003) and Sargent et al. (2006), among others. Additionally, we do not 
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observe crystalline grains with sizes larger than 1.5 fim, in agreement with Bouwman et al. (2006), 
who suggest that this fact may be due to agglomeration of the (less abundant) crystalline silicates 
with the bulk of amorphous grains, producing mixed crystalline and amorphous aggregates, rather 
than large, fully crystalline grains. The lack of correlation between age and crystallinity suggests 
that crystallinization may occur already at the very early stages of star formation and disk evolution 
(Dullemond et al. 2006). 

We have IRS spectra for two "transition objects" with accretion rates below 10~ 12 M /yr 
(consistent with no accretion, confirmed via both the lack of an U band excess and the narrow 
Ha profiles at R~34,000 resolution, Paper IV), named 73-758 and 14-11, for a transitional disk 
with ongoing accretion (24-515), and for nearly-transitional disk (a disk with nearly photospheric 
near-IR colors), 13-1250. The silicate feature is very weak in 14-11, and absent or extremely faint 
in 24-515, so only 73-758 could be fitted. Although we cannot fit the silicate feature of some 
"transition objects" , its weakness is consistent with the general trend showing that stars with lower 
accretion rates tend to have weaker silicate features, and supports the idea of dust evolution and 
gas evolution happening somehow in parallel in accretion disks. 



4.2. Age, Inner Disk Evolution, and the Silicate Feature 

In order to complete our understanding of the differences in silicate emission due to age and 
spectral type in these stars, we have analyzed the IRS spectra grouping them in classes according 
to similar ages and similar spectral types. Figures [18] and [19] display the younger population (0-2 
Myr, all but 91-155 associated with the young globule population, see Paper II), the average age 
in TV 37 (2-6 Myr), and the older group members (over 6 Myr, including the only accretor in NGC 
7160, 01-580, aged 9 Myr). Figures [201 and [211 display the spectra in groups with similar spectral 
types (G9-K2, K4-K5, K6-K7, and M0-M2.5), containing (except the first one) about the same 
number of stars in each group. 

As for the age classification, we observe a large variety in the silicate emission morphology and 
strength. Individual ages of young stars must be handled with care, as some errors may occur due 
to the presence of variability and the problems defining the birthline (affecting mostly G stars, see 
Hartmann 2003). Nevertheless, studies including variability effects (Burningham et al. 2005) have 
revealed that ages derived from the V vs. V-I diagram are not much affected by variability, so that 
the main differences observed in these color-magnitude plots are caused by real age differences. 
Therefore, as discussed in Paper II, our ages are likely to be accurate up to ~1 Myr, except for 
the G stars (the star 82-272 in our sample), for which the uncertainties are larger. As noted in 
Section 14. 1[ we observe that half of the stars with atmospheric silicates dominated by amorphous 
small grains are among the older group with ages over 6 Myr. This is the case of the NGC 7160 
star 01-580, and of the 7 Myr-old 23-162. There are only 2 younger stars (among 16 with ages <6 
Myr), 12-2113 (a SB1, aged 1 Myr), and 13-157 (aged 2 Myr), which contain mostly unprocessed 
small amorphous silicates (with average sizes <2^m) similar to 01-580 and 23-162. Despite the 
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low-number statistics, the proportion of stars with pristine silicate features seems higher for older 
stars, which could be an additional indication of dust settling removing the large grains from the 
upper layers of the disks with time, in agreement with the anti-correlation between size and age 
described in Section 14.11 

The classification according to the spectral types shows also interesting results that escaped 
our analysis from the silicate fitting, given that 7 of the low-mass spectra with very weak silicate 
features could not be properly fitted. Observing Figures [20] and EH we can see that, despite the 
large variety of silicate features present in all the spectral type groups, all the spectra with absent 
or nearly absent silicate features belong to the later spectral types, being 2 out of the 7 spectra with 
K6-K7 types, and 4 of the 5 spectra with M0-M2.5 types. Among the 8 spectra with spectral types 
earlier than K6, we find only one example of small but certainly well-detectable silicate feature, 
21-2006. The peak of the flux over the normalized continuum of the K5 star 21-2006 is 1.7±0.1, 
larger the peak fluxes of stars like 12-1209, 11-2322, 14-11, 11-2131, 24-515, and 21-33, ranging 
from 1.2 to 1.5±0.1. A KS test comparing the normalized peak flux of stars with spectral types 
earlier than K6, and K6 and later, gives a probability of 0.03 (3%) that both samples come from 
a similar distribution. Taking into account that we have 8 spectra earlier than K6, and 11 spectra 
with types K6 or later, and that they were drawn randomly from our member sample in Cep OB2, 
this difference has a statistical significance. About half of the late-type stars have very weak silicate 
features, compared to roughly 10% of the earlier-type stars, or said in a different way, nearly 90% 
of the stars with very weak or absent silicate feature in our sample are K6 or later. 

If we examine the sample of TTS studied by Honda et al. (2006), we find the same behavior: 
there is only 1 star with very weak silicate feature among their 8 stars with spectral types G2 to 
K5, but there are 6 among their 22 stars with spectral types K6 to M4 (all of them M-type, except 
for one K6 star), so ~85% of their stars with weak silicates are K6 or later. Although our sample 
seems richer in stars with weak silicate features (about 50% of the later spectral types, compared 
to 30% in Honda et al. 2006), the trend is the same: Stars with very weak or undetectable silicate 
feature are 3-5 times more frequent among the spectral types later than K6. Honda et al. (2006) 
relates the lack of silicate to strong accretion, which is obviously not the case in our sample (we 
moreover observe the opposite, as stars with stronger accretion rates tend to have stronger silicate 
features). Nevertheless, all but one of the stars in their sample have Ha normal for low-mass young 
stars accreting at typical rates, taking into account that M stars have larger Ha EW than K stars 
for similar accretion rates (White & Basri 2003). The star DN Tau (that they mention as having 
a large accretion rate causing the reduced silicate emission) was found to have an accretion rate of 
3.5x lO _9 M0/yr (Gullbring et al. 1998), so the presence of systematically higher accretion rates is 
not clear. Moreover, Green et al. (2006) has shown that even FU Ori objects have strong silicate 
features, despite their extremely high accretion rates. Most of the stars in the sample of Honda 
et al. (2006) are Taurus stars, younger than ours, although they include a few from the TW Hya 
association (10 Myr age). If inner disk evolution increases with time, as suggested in Paper III, 
we would then expect that the number of stars with weak silicates increases with time. But even 
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considering our sample, it is not easy at this point to derive a correlation between age and the 
proportion of stars with no silicates due to the few observations of disks aged 10 Myr. In the 
sample of stars investigated by Kessler-Silacci et al. (2005), which contains mostly younger stars as 
well, the M-type stars tend to have smaller peak fluxes than the higher-mass ones, although their 
sample has very few stars with weak silicate features, and some spectral types are largely uncertain. 

During some time, there was the open question of whether late M stars and brown dwarfs 
(BD) would show any silicate emission, precisely because of these requirements for having warm 
silicates. Nevertheless, the work of Apai et al. (2004, 2005) and Mohanty et al. (2004), among 
others, revealed that BD can have strong silicate features, including strongly processed grains, and 
there are no doubts that early M-type stars as those in our sample can have strong silicate emission 
(Bouwman et al. 2001, 2006; Sargent et al. 2006; among others). But most of the systems studied 
contain younger stars (typically, 1-2 Myr), and this may explain part of the difference, taking into 
account the changes in stellar luminosity and the evolution of the inner disk with age. 

There are several factors that could be responsible for the very weak silicate features in M- 
type stars, and one of them could be the region where the silicate feature is produced, the warm 
disk atmosphere. Due to the lower effective temperature of M-type stars, the silicate feature is 
produced in a region much closer to the star (and thus smaller) in M stars than in early K and 
G stars. In order to achieve a temperature around 300 K, required for the silicate emission at 10 
^m, the small grains in the atmosphere need to be at distances closer than 0.7 to 1.5 AU for an 
M0 star, compared to 1.2 to 3 AU for a K5 star, approximately, and the area of the disk capable 
of producing silicate emission is therefore about 4-5 times larger in a K5 star than in a M0. As the 
luminosity of the star decreases with age, the warm silicate region moves inwards. The change of 
luminosity is remarkable for M-type stars, for which it may drop by an order of magnitude from the 
age of 1 to 10 Myr. For a typical M2 star, the warm disk region would change from being located 
at ~0.4 to 1.2 AU, to being reduced to ~0.1 to 0.4 AU. Moreover, we observed that important 
inner disk evolution has occurred by the ages of ~ 4 Myr (Paper III), which can be inferred from 
the weaker near-IR emission, suggesting significant grain growth and/or dust settling. Therefore, 
we may not only be probing a different area of the disk, but a more evolved one as well, if dust 
evolution (dust coagulation/settling/opacity changes/gap opening) is faster and more efficient in 
the innermost disk. Although we do not observe any significant difference between the inner disk 
emission, including the presence of inner gaps, for the different spectral types from late G to early 
M (Paper IV), grain growth and/or dust settling would affect the silicate emission of M stars more 
than for K stars. K and G stars could have part of their inner disks depleted of small grains (as it 
may happen to the K4.5 star 13-1250, which has a very reduced near-IR excess that makes it very 
close to a "transition object"), but the higher stellar luminosity will be still enough to produce 
the temperatures required for silicate emission at larger distances from the disk, unlike in M-type 
stars. Therefore, the silicate feature could be a good instrument to probe the evolutionary stage of 
the innermost disk in late type stars, even before other effects like SED flattening in the near-IR 
become evident. 
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Although our study in Papers III and IV did not reveal any correlation between the presence 
of inner gaps and the spectral type, we noted in Paper IV that other works have suggested that 
near-IR gaps (detectable at IRAC wavelengths) tend to appear with a higher frequency around 
M-type stars (McCabe et al. 2006), which could be a proof of photoevaporation removing the 
inner disk (the gaps around earlier-type stars would be larger, and not detectable with IRAC). Cep 
OB2 is again the older region where a statistically significant study of transitional disks has been 
made. Therefore, the fact that at 4 Myr we do find inner gaps in K stars similar to the gaps found 
for M star by McCabe et al. 2006 at 1-2 Myr may as well indicate that the evolution of M-type 
stars differs in some aspects from the evolution of later spectral types. Nevertheless, since we find 
that half of the transition disks are still accreting, we need probably other effects instead of or in 
addition to photoevaporation, in order to produce the inner gaps. The idea of M stars suffering 
evolutionary processes different from the processes affecting their higher mass counterparts could 
result in the differences observed here in the strength and shape of the silicate emission. Lower 
silicate emission might as well suggest less turbulence in the lower-mass disks, and maybe different 
disk structure and evolution. If disks around M-type stars differ in some aspects from more massive 
T Tauri disks, it could explain the difficulties for M-type stars to form giant planets suggested by 
the models of Laughlin et al. (2004) and the lack of giant planets around the less massive stars in 
the observations of Endl et al. (2006). Despite the historical considerations that K and early M 
stars evolve similarly, recent studies suggest differences in disk evolution around M stars, detectable 
as lower disks fractions and higher percentages of "flattened" disks at a given age, compared to 
K-type stars (Lada et al. 2006). In addition, accretion mechanisms operating in M-type stars and 
very low-mass objects may also need to be reconsidered, as the depth of the active layer of the disk 
and the presence of dead zones may be very different in case of the very low-mass disks around 
late- type objects (Hartmann 2006). Further observations of very low-mass stars at different ages 
are strongly required in order to investigate all these hypotheses. 



4.3. Transitional and debris disks around intermediate- and high-mass stars 

The IRS data on high- and intermediate-mass stars in Tr 37 reveal that two of these stars 
(CCDM+5734Ae/w and KUN-196) are more likely non-accreting transitional protoplanetary disks 
with inner gaps, as the silicate feature at 10/im is present. The 10/um feature suggests that, despite 
of the inner gaps, some small grains are present at distances of the order of 10-50 AU for the B3+B5 
star CCDM+5734Ae/w, and from -5 to 15 AU for the B9 star KUN-196. KUN-314s may be either 
an evolved Herbig Ae star (with a disk flatter than a typical Herbig Ae star), or a classical Ae 
star. The rest of the stars (MVA-1312, MVA-447, MVA-468, BD+752356) do not show any silicate 
emission, so they are consistent with typical "debris disks" as stated in Paper III. In NGC 7160 
we do not find any Herbig Ae/Be star. The A7 star DG-481 shows a very clear silicate feature, 
produced mainly by small amorphous grains that would be located at distances between ~2 and ~6 
AU, so the system is probably a non-accreting transitional protoplanetary disk. The rest of stars 
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(DG-39, DG-682 and DG-9l4j) do not show any evidence of silicate emission (Figure 0|), being 
consistent with "debris disks" (Sylvester & Mannings 2000). 

Analyzing the high- and intermediate- mass disks, we must mention the difference in proto- 
planetary disk frequencies compared to low-mass stars, given that only 3-4 stars in Tr 37 (out of 
58 members) and 1 star in NGC 7160 (out of 66 members) present protoplanetary disks (with or 
without inner gaps). The lack of classical Herbig Ae/Be stars in the older cluster and the fact 
that the total number of stars with excesses is double in Tr 37 compared to NGC 7160, due to 
the presence of more potential protoplanetary disks with large gaps in Tr 37, seems to point to 
disk evolution among these higher-mass stars (given that they contain 58 and 66 ABF members, 
respectively). 

The fraction of potential debris disks is similar in both clusters, which would be also the 
expected result if the signatures we observe correspond to the disks left over and stirred by recently- 
formed planetary systems (Kenyon & Bromley 2005). In Table [6] we list the excesses over the 
photospheric levels of all the high- and intermediate-mass stars. In agreement with Kenyon & 
Bromley (2005), the observed excess luminosity over the photosphere of the potential debris disks 
is larger than the luminosity seen in older systems. This is consistent with the predictions of the 
models of debris disk evolution, which include large bodies (>1000 km planets/planetesimals) in 
the inner disk stirring the dust and replenishing the debris disk via collision events. The excess 
flux is lower than the extremely bright HR 4796A disk, which is brighter than the predictions of 
typical young debris disk models. 

According to Rieke et al. (2005), we would expect an excess at 24 /im of the order of 1.25 to 5 
at ages under 24 Myr, although the fraction of stars with a strong excess (~5) would be larger (11% 
compared to 5%) for objects under 10 Myr age. The decay of the luminosity would be exponential, 
with a decay time of the order of hundreds of years. According to Kenyon & Bromley (2005), 
the predicted typical excesses over the photosphere at 24//m for an A star, after planet formation 
in the innermost disk (<3 AU), would be of the order of log(F24/F24 P )~1.25-0.8 at 1 Myr after 
the formation of the planets, 0.8-0.5 after 4 Myr, and 0.2-0.6 after 12 Myr. Nevertheless, if the 
planets were to be formed at 3-20 AU, the excess could be of the order of ~1.8 after 1 Myr, to 
decay to ~1 after 10 Myr, and even 1.2 after 10 Myr if the planets are formed at 30-150 AU. These 
predictions are consistent with our observations, so all the disks but those of MVA-426 and KUN- 
314S are consistent with debris disks luminosities, as stated in Paper III. Nevertheless, the presence 
of silicates in CCDM+5734A, KUN-196 and DG-481 suggest that they are rather protoplanetary 
disks with large gaps. Although our sample is limited by the small number of stars, debris disks in 
Tr 37 tend to be brighter than their counterparts in NGC 7160 (as we would expect from debris 
disk evolution), but here we should take into account that the stars in NGC 7160 have later spectral 
types. 



2 The excess in DG-62 is probably not significant, since it is consistent with photospheric levels within the errors, 
so we cannot confirm the detection of a disk for this star. 
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However, it would be necessary to examine these disks at longer wavelengths (mm and submm) 
in order to determine if their opacity is mostly due to second-generation dust from collisions, and 
if they are optically thin, as the weak excess seem to suggest in all or part of the cases (see Table 
[6]). Longer wavelength observations are specially important, given that in our low-mass sample we 
have observed several examples of stars with optically thick protoplanetary disks at 24 //m and 
extremely weak near-IR excess and silicate features. The higher temperatures of the central stars 
and the larger areas subtended by the inner disks here would require much larger inner gaps in 
order to avoid silicate emission, making it more plausible that they are at least partially optically 
thin. In that case, the stars detected in Tr 37 and NGC 7160 would be among the youngest debris 
disks detected, and some of the few ones with well determined ages (derived from the average ages 
of the low-mass clusters members), showing that debris disks are present already at ages as young 
as 4 Myr. The number of debris disks observed at 4 and 12 Myr is roughly similar (4 plus maybe 1 
versus 3 plus 1 uncertain, respectively). The emission over the photosphere of the presumed debris 
disks is variable in both clusters, as we would expect if initial and/or environmental conditions (and 
maybe the characteristics of a planetary system that it could have formed there) play an important 
role in disk evolution, in addition to differences in age and mass (Kenyon & Bromley 2005; Rieke 
et al. 2005). 



5. Conclusions 

We presented IRS spectra of all the 13 high- and intermediate-mass stars with IR excess, 
as well as a total of 20 low-mass stars selected among the stars with disks in the clusters Tr 37 
and NGC 7160. Our sample contains stars that are on average older than the members of the 
best known regions, so we can examine with more detail the effects of disk evolution in the shape 
and strength of the silicate emission. The IRS data reveal a large variety of spectra and silicate 
features in both clusters, suggesting different stages in dust evolution, coagulation and presence of 
small-grains in the warm disk atmosphere, and probable dependence on initial conditions and/or 
environment (Przygodda et al. 2003; Kessler-Silacci et al. 2005; Dullemond & Dominik 2005). We 
do not find any significant differences in grain size nor composition between the low-mass (G-M) 
stars and the only early-type star with strong 10/um feature, revealing that disks, silicate grains, and 
dust processing can be similar in systems with central stars of very different masses. Nevertheless, 
despite of the large variety of cases, there are significant trends that relate the disk structure, 
turbulence, the amount of settling, and the spectral type to the strength and characteristics of the 
silicate feature. 

We observe clear evidence of the dependence of silicate emission with spectral type. Objects 
with very weak silicate features are 3-5 times more frequent among the M-type stars than among 
K-type. The differences may be larger if we consider early K and G-type stars, as we do not 
observe any G or early K stars with very weak silicate features, although the total numbers are 
quite similar. One possibility to explain this could be inner disk evolution, given that the silicate 
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feature is produced closer to the star and in a smaller region in M-type star than in K and earlier 
types. Therefore, we would be probing different disk regions in K and M stars, which may have 
suffered differential disk evolution (grain growth/dust settling). The probed regions would change 
with time as the stellar luminosity decreases. In that case, the silicate feature may be a very good 
indicator of the evolutionary stage of the inner disk, even before than the effects of the grain growth 
and the sedimentation produce detectable disk flattening. It could also indicate that evolutionary 
processes and timescales are somehow different for the later spectral types, as some authors have 
suggested (based on disk fractions and accretion; Lada et al. 2006; Hartmann 2006), and this could 
have effects in their capabilities to form planetary systems. 

We find an apparently counter-intuitive correlation that large average grain sizes (~5-6 //m) 
are found only among the youngest stars (aged 0-2 Myr) and the stars in our sample that show very 
pristine silicate features (containing only small amorphous grains) tend to be among the oldest ones 
(>6 Myr). This trend is consistent with dust settling increasing with time, as large grains would 
drift faster towards the middle plane, leaving only very small and unprocessed grains in the disk 
atmosphere at the later stages of disk evolution. Additionally, the fact that we observe already large 
grains at very early ages suggests that considerable grain growth occurs already within the very 
first years of disk formation and evolution, as suggested by Dullemond et al. (2006). Therefore, 
one must be cautious before assigning an earlier evolutionary stage to objects with silicate features 
due to smaller amorphous grains. 

We observe that stars with low accretion rates have in general smaller grains than the strong 
accretors, which could suggest that significant dust settling has occurred in the more evolved, slow 
accreting, stars, and/or that the presence of turbulence in disks would play an important role 
replenishing the disk atmosphere from silicate grains that may otherwise settle into the middle 
plane with time. This could be an additional evidence of dust evolution occurring (or requiring) 
parallel to gas and accretion evolution (as we mentioned in Paper IV), as it may be required that 
gas accretion decreases in order to allow the dust settling to take place. 

In the case of transitional disks with inner gaps, we find very weak features among the 3 
objects examined here, consistent with the trend of weak silicate emission related to stars with 
very small or zero accretion. We find that one of the transition objects, the Ml. 5 star 14-11, is 
very similar to CoKu Tau/4, having no accretion and emission at A >8 fim consistent with a wall 
or rim surrounding an inner gap approximately 10 AU in size. 

Regarding crystalline silicates, we confirm the observations of Meeus et al. (2003), Bouwman et 
al. (2006), and Honda et al. (2006), finding that there is no clear trend in the degree of crystallinity 
versus spectral type, age, accretion rate nor SED slope. We find that crystalline grains tend to be 
smaller than amorphous ones, and that they are about the same size at any age or evolutionary 
stage, maybe because of the difficulties to form entirely crystalline large grains (Bouwman et al. 
2006). Crystalline grains may be therefore part of aggregates with the more abundant amorphous 
silicates. The degree of crystallinity is in general low, being always well below 20%. The two 
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objects with higher (>15%) crystallinity fractions are two spectroscopic binaries, which would be 
in agreement with Meeus et al. (2003), who reported that binaries tend to have more crystalline 
disks, probably because their disks are more turbulent. 

Among the high-mass stars, we can confirm a total of 7 potential debris disks objects with 
average ages 4 and 12 Myr. These potential debris disk lack the silicate emission and show mid-IR 
excess comparable to other typical debris disks, and consistent with the excess predicted a few Myr 
after planet formation has occurred (Kenyon & Bromley 2005; Rieke et al. 2005). Further research 
at longer wavelengths (mm and submm) would be required to check whether these disks are debris 
disks, or they are in a more evolved phase within the transitional stage from Herbig AeBe disks to 
Vega-like systems. 

The correlations between accretion rates/turbulence, spectral type and the silicate feature 
observed here can fit within the same picture: With age, both the gaseous and the dusty components 
of the disk evolve in a parallel way. The decrease of the accretion rate (consistent with viscous disk 
evolution) is accompanied by the settling of dust grains to the mid-plane, and the lower turbulence 
levels are not enough to replenish the disk atmosphere with more small grains. The settling is 
stronger in the innermost disk, and would extend outward with time, until strong dust coagulation 
results in the opening of opacity gaps and eventually, maybe in planet formation. 

We want to acknowledge B. Merin and C. Dullemond for interesting discussion. We also want 
to thank the referee for the comments and suggestions to clarify this paper. 
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Fig. 1. — Spectral Energy Distribution of low-mass stars in Cep OB2. Optical (UVRI), 2MASS 
(JHK), IRAC (3.6-8.0 /xm) and MIPS (24/zm) data is represented by dots. The IRS spectrum is 
displayed with its error. A photosphere from Kenyon & Hartmann (1995) is displayed for each 
spectral type. The flux in log(A Fa) is displayed in units of erg cm~ 2 s _1 . 



-26- 




Fig. 2. — Spectral Energy Distribution of low-mass stars in Cep OB2 (continued). Caption as in 
Figure CD 
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Fig. 3. — Spectral Energy Distribution of intermediate- and high-mass stars in Tr 37. Optical 
(UB), 2MASS (JHK), IRAC (3.6-8.0 fan) and MIPS (24pm) data is represented by dots. The IRS 
spectrum is displayed with its error. A photosphere from Kenyon & Hartmann (1995) is displayed 
for each spectral type. The flux in log(A Fa) is displayed in units of erg cm -2 s" 1 . 



-28- 




Fig. 4. — Spectral Energy Distribution of intermediate- and high-mass stars in NGC 7160. Optical 
(UBV), 2MASS (JHK), IRAC (3.6-8.0 /im) and MIPS (24pm) data is represented by dots. The IRS 
spectrum is displayed with its error. A photosphere from Kenyon & Hartmann (1995) is displayed 
for each spectral type. The flux in log(A F>) is displayed in units of erg cm~ 2 s _1 . 
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Fig. 5. — Fit for the Silicate feature at 10/um for low-mass stars in Cep OB2. The observed 
spectrum is displayed with its errors, and the best fit is drawn as a thick line. The continuum level 
is displayed as a dotted line; the emission due to amorphous grains is displayed as a dashed line; 
the emission due to crystalline grains is represented by a dotted-dashed line; the PAH contribution 
is a long-dashed line. Vertical lines are displayed to mark the locations of the peaks of enstatite 
(9.2 /mi), amorphous silicate of olivine and forsterite stoichiometry (9.7 and 10.2 /mi, respectively), 
and crystalline forsterite (11.3 /mi). 
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Fig. 6. — Fit for the Silicate feature at 10/im for low-mass stars in Cep OB2 (continued). See 
caption in Figure [5j 
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Fig. 7. — Fit for the Silicate feature at 10//m for low-mass stars in Cep OB2 (continued). See 
caption in Figure 
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Fig. 8. — Fit for the Silicate feature at 10/um for intermediate- and high-mass stars in Tr 37. The 
observed spectrum is displayed with its errors, and the best fit is drawn as a thick grey line. The 
continuum level is displayed as a dotted line; the emission due to amorphous grains is displayed 
as a dashed line; the emission due to crystalline grains is represented by a dotted-dashed line; the 
PAH contribution is a long-dashed line. Vertical lines are displayed to mark the locations of the 
peaks of enstatite (9.2 /jm), amorphous silicate of olivine and forsterite stoichiometry (9.7 and 10.2 
fim, respectively), and crystalline forsterite (11.3 jim). In the case of the Herbig B7 star MVA-426, 
which lacks any silicate feature, we mark the strong PAH features at 6.2, 7.7, 8.6, 11.3 and 12.7 
fj,m. 
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Fig. 9. — Fit for the Silicate feature at 10//m for intermediate- and high-mass stars in NGC7160. 
The observed spectrum is displayed with its errors, and the best fit is drawn as a thick grey line. 
The continuum level is displayed as a dotted line; the emission due to amorphous grains is displayed 
as a dashed line; the emission due to crystalline grains is represented by a dotted-dashed line; the 
PAH contribution is a long-dashed line. 
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Fig. 10.— SEDs of the 3 "transition disks" 73-758, 14-11, and 24-515, and the very evolved disk 
13-1250. The optical (UVRI), JHK (from 2MASS), IRAC and MIPS (24 pm) photometry points 
are displayed together with the IRS spectrum. The stellar photosphere and the median CTTS 
SED from Taurus (Hartmann et al. 2005) are shown for comparison as the dotted-dashed line and 
the dashed line, respectively. The Ha line (observed with Hectochelle, R~ 34,000) is shown in the 
figure inset. Note that 73-758 and 14-11 do not show any velocity wings, and that the velocity 
wings of 24-515 are only marginal. 
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Fig. 12. — Ratios of the normalized flux in the 10 /jm region for the stars showing silicate features 
in TV 37 and NGC 7160, compared to the peak of the normalized flux and to two stellar properties 
(spectral type and age). The typical error bar is shown. The A7 star DG-481 is marked with a 
large square to distinguish it from the rest of the sample, containing only T Tauri stars. Note that 
DG-481 is displayed with a shift in the spectral type plot. The trend for amorphous grains with 
sizes increasing from 1 to 10/xm is shown in the top left panel as a dotted line. The data points of 
Bouwman et al. (2006) are represented as well in the panel, marked by crosses. In the top right 
figure, the trends of the data of Bouwman et al. (2006) for the F11.3/F9.8 and Fg.e/Fg.s are shown 
as dotted and dashed lines, respectively. 
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Fig. 13. — Peak of the normalized flux in the 10 /im region for the stars showing silicate features 
in Tr 37 and NGC 7160, versus different disk and stellar properties (accretion rate and IR colors). 
The A7 star DG-481 is marked with a box. Typical error bars are displayed in each diagram. The 
correlation between the flux at the normalized peak and the [24]-[31] color is consistent with the 
data of Bouwman et al. (2006) for low-mass stars, shown as crosses in the bottom right panel. The 
trend showing that earlier spectral type stars tend to have higher accretion rates is shown in the 
top right panel for comparison with the trends relating accretion rates to other silicate and stellar 
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Fig. 14. — Grain size versus some stellar and disk properties (spectral type, age, and accretion 
rate). The A7 star DG-481 is marked with a large square (note that it is displaced in the spectral 
type figure for clarity), the rest of stars are TTS. The fact that larger (amorphous) grains are found 
preferentially in young stars suggest dust settling with age. Most of the large grains are found in 
stars with the higher accretion rates, as it would be expected if the disk atmosphere is replenished 
by stirring the grains that may be settling in more turbulent disks. The described behavior seems to 
affect only the amorphous grains (which constitute most of the mass), the behavior of the (smaller) 
crystalline grains is not found to correlate significantly with other properties than binarity. The 
correlation between grain size and normalized flux strength (F pea k) is also displayed. 
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Fig. 15. — Grain size versus the SED slope at different radii. The A7 star DG-481 is marked with 
a large square, the rest of stars are TTS. 
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Fig. 16. — Total mass from the silicate fit, and mass in amorphous grains, versus different disk and 
stellar properties (spectral type, age, accretion rate, and near-IR color). There is a weak correlation 
between the silicate mass and the H-K color (related to the flaring of the inner disk/inner wall), 
suggesting that the amount of warm silicates depends on the irradiation (therefore, temperature) 
of the inner disk. No correlation is observed between mass and spectral type, as we would expect if 
the amount of silicates depends rather on disk geometry and properties than on stellar properties. 
The correlation between the mass in silicates and the accretion rate could be related to the higher 
turbulence of the stronger accretors, but could be also due to the differences in spectral type. The 
A7 star DG-481 is marked by a large square. Note that its spectral type is shifted for better viewing 
of the low-mass stars. 
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Fig. 17. — Fraction of crystalline silicates versus different characteristics of the disks and stars. The 
fraction of crystalline silicates does not seem to be correlated with any stellar property (spectral 
type, SED shape, mass accretion rate, age), although it is remarkable that the two stars with larger 
crystalline fraction are spectroscopic binaries. The A7 star DG-481 is marked by a large square. 
Note that its spectral type is shifted for better viewing of the low-mass stars. 
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Fig. 18. — IRS spectra for stars with ages 0-2 Myr (left) and 2-6 Myr (right). The location of 
the main silicate peaks (9.7 and 18 fim peaks of the amorphous silicates of olivine stoichiometry, 
10.2 and 23.5 //m peaks of the amorphous silicates of forsterite stoichiometry, 11.3 /xm crystalline 
forsterite, and the 28.5 and 33.5 jim. crystalline bands) is marked by dashed lines; the SED from 
our photometry studies are depicted as well. Spectra are shifted in flux to allow the visualization 
of them all, as absolute fluxes do not matter for our discussion. 
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Fig. 19. — IRS spectra for stars with ages over 6 Myr. Note that the age for the G9 star 82-272 
is uncertain. As in Figure [TE[ we have shifted the spectra for display and we mark the location of 
the main silicate peaks by dashed lines. 
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Fig. 20. — IRS spectra for stars with spectral types G9-K2 (left) and K4-K5 (right). As in Figure 
[T8| we have shifted the spectra for display and we mark the location of the main silicate peaks 
by dashed lines. Among the earlier spectral types, the lack of silicate emission seems to be highly 
infrequent. 
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Fig. 21. — IRS spectra for stars with spectral types K6-K7 (left) and M0-M2.5 (right). As in 
Figure [THl we have shifted the spectra for display and we mark the location of the main silicate 
peaks by dashed lines. Stars with later spectral types seem to display a lack of silicate emission 
more frequently (see discussion in text). 



Table 1. Summary of Disk and Stellar Properties of the Observed Stars 



Cluster 


ID 


RA(2000) 


DEC(2000) 


Sp.Type 


Age (Myr) 


M (10- 8 M Q /yr) 


Class 


Comments 


Tr 37 


73-758 


21:35:08.39 


57:36:03.2 


K6.5 


1.8 





TO 




Tr 37 


11-2146 


21:36:57.69 


57:27:33.3 


K6 


0.9 


16.2-13.2 


II 




Tr 37 


11-1209 


21:36:58.55 


57:23:26.1 


K6 


1.0 


<0.1 


II 


10/im feature marginal 


Tr 37 


11-2322 


21:37:01.96 


57:28:22.3 


Ml 


0.8 


0.80 


II 


10/im feature marginal 


Tr 37 


11-2037 


21:37:07.07 


57:27:01.1 


K4.5 


2.5 


0.97-2.5 


II 


28/im Forsterite 


Tr 37 


14-11 


21:37:10.30 


57:30:19.1 


Ml. 5 


0.7 





TO 


~10 AU gap, BB~70 K, 10/im feature marginal 


Tr 37 


11-2131 


21:37:12.19 


57:27:26.6 


K6.5 


2.3 


1.1 


II 


28/im Forsterite, 10pm feature marginal 


Tr 37 


11-2031 


21:37:15.95 


57:26:59.5 


K2 


2.5 


1.6 


II 




Tr 37 


14-183 


21:37:38.53 


57:31:41.7 


K6 


0.9 


<0.1 


II,SB1: 




Tr 37 


82-272 


21:38:03.53 


57:41:35.4 


G9 


10.5 


23.9 


II+II,SB2 


18/im Silica, 23-28/im Forsterite 


Tr 37 


12-2113 


21:38:27.47 


57:27:21.3 


K6 


1.1 


6 


II,SB1 




Tr 37 


13-157 


21:38:28.12 


57:30:47.1 


K5.5 


2.4 


1.6 


II 




Tr 37 


91-155 


21:38:34.75 


57:41:28.0 


M2.5 


1.7 




II 


i 
i 


Tr 37 


54-1547 


21:38:44.48 


57:18:09.6 


K5.5 


5.7 


0.212 


II 




Tr 37 


13-1250 


21:39:12.17 


57:36:16.9 


K4.5 


3.3 


0.10 


TO/II 


Small excess in near-IR, strong settling/growth 


Tr 37 


24-515 


21:39:34.10 


57:33:32.1 


M0. 5 


2.8 


<0.1 


TO 


10pm feature marginal 


Tr 37 


21-33 


21:39:35.64 


57:18:22.3 


M0 


4.0 


<0.1 


II 


10pm feature marginal 


Tr 37 


21-2006 


21:40:13.91 


57:28:48.5 


K5 


6.1 


0.09 


II 




Tr 37 


23-162 


21:40:44.52 


57:31:31.9 


K7 


6.6 


<0.1 


II 




NGC 7160 


01-580 


21:53:37.07 


62:28:47.2 


K4.5 


8.7 




II 




Tr 37 


MVA-1312 


21:35:56.82 


57:20:53.1 


B4 


4( ) 





D 


rirvn ATT Tl Tl 1 /\r\ f f 

~300 AU gap, BB~100 K 


Tr 37 


KUN-314s 


21:36:14.45 


57:21:32.5 


A: 


4(*) 




TO :/cBe: 


Featureless 


Tr 37 


CCDM+5734A 


21:37:40.93 


57:33:37.5 


B3+B5 


4(*) 





TO, VB 


18/im Silica, Gap with small grains within 10-50 


Tr 37 


MVA-468 


21:38:08.18 


57:31:26.7 


B7 


4 (*) 





D 




Tr 37 


MVA-426 


21:38:08.71 


57:26:48.0 


B7 


4(*) 




HAeBe 


PAH 


Tr 37 


MVA-447 


21:38:26.64 


57:28:40.8 


F0 


4(*) 





D 




Tr 37 


BD+572356 


21:38:30.31 


57:46:26.5 


A4 


4 (*) 





D 




Tr 37 


KUN-196 


21:40:15.14 


57:40:51.2 


B9 


4(*) 





TO: 


10/1 feature marginal, small grains at ~5-15 AU 


NGC 7160 


DG-481 


21:52:21.10 


62:45:03.7 


A7 


12(*) 





TO 


10/1 feature strong, small grains at ~2-6 AU 


NGC 7160 


DG-39 


21:53:27.73 


62:35:19.0 


AO 


12(*) 





D 





Table 1 — Continued 



Cluster 


ID 


RA(2000) 


DEC(2000) 


Sp.Type 


Age (Myr) 


M (lO- 8 M /yr) 


Class 


Comments 


NGC 7160 


DG-682 


21:53:45.03 


62:36:55.2 


A2 


12(*) 





D 




NGC 7160 


DG-62 


21:54:30.23 


62:31:16.1 


F5 


12(*) 





D: 


Excess not significant 


NGC 7160 


DG-912 


21:55:47.50 


62:35:43.7 


F5.5 


12(*) 





D 





Note. — Summary of the properties of the stars observed with IRS. For the low-mass stars, ages are derived from the V vs. 
V-I diagram, using the isochrones by Siess et al. 2000 (see Sicilia-Aguilar et al. 2005 for a detailed description of spectral typing 
and age/mass calculation). Accretion rates are based on U band emission (see Sicilia-Aguilar et al. 2005, 2006a, 2006b) or are 
upper limits from high-resolution Ha spectra. We denote with M=0 those stars with no U excess and Ha profiles ( M<10~ 12 
M©/y r > Paper IV). The "Class" can be class II (II), class III (III), transition object (TO), Herbig Ae/Be star (HAeBe), classical 
Be star (cBe) or debris disk (D). We also include notes about binarity (SB1= single-lined spectroscopic binary; SB2=double- lined 
spectroscopic binary; VB=visual binary). The "Comments" column include notes about the IRS spectra, among others. (*) For 
the high- and intermediate-mass stars, no individual ages nor masses are estimated, and we do not have any evidence of accretion. 
Uncertain values are denoted by ":". 
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Table 2. Scale Factors for the Combination of Orders 



ID 



/ o- ( Oo 


IN O 


scaling 


11-2146 


SL 


x 1.20 


11-1209 


No 


scaling 


11-2232 


No 


scaling 


11-2037 


No 


scaling, 


14-11 


No 


scaling 


11-2131 


No 


scaling 


11-2031 


No 


scaling 


14-183 


No 


scaling, 


82-272 


No 


scaling 


12-2113 


No 


scaling 


13-157 


No 


scaling 


91-155 


No 


scaling 


54-1547 


No 


scaling, 


13-1250 


No 


scaling 


24-515 


No 


scaling 


21-33 


No 


scaling 


21-2006 


No 


scaling, 


23-162 


No 


scaling, 


01-580 


No 


scaling 


MVA1312 


No 


scaling, 


KUN-314S 


No 


scaling 


CCDM+5734Ae 


SL 


x 1.1 


MVA-468 


No 


scaling 


MVA-426 


No 


scaling 


MVA-447 


No 


scaling 


BD+572356 


No 


scaling 


KUN-196 


SL 


x 6.50, 



used LL bonus order from 19.40 to 21.00 fim 



used LL bonus order from 19.40 to 21.00 /im 



from 19.40 to 21.00 £tm (LL) 



LL bonus order from 19.40 to 21.00 //m 



DG-481 

DG-39 

DG-682 

DG-62 

DG-912 



No scaling, used bonus order from 7.53 to 8.10 (im (SL), from 19.40 to 21.00 (im (LL) 
SL x 4.74, used bonus order from 7.53 to 8.10 ^m (SL), from 19.40 to 21.00 /im (LL) 
SL x 2.00, used bonus order from 7.53 to 8.10 ^m (SL), from 19.40 to 21.00 ^m (LL) 
SL x 1.717 

SL2 x 1.29, SL1 x 1.79, LL1 x 2.48 



Note. — Scale factors in the data reduccion process. Where indicated, the SL or LL orders were multiplied 
by a certain scaling factor in order to match the other orders. 



Table 3. Dust species used to model the silicate feature 



Species 



State 



Shape 



Size (fim) 



Reference 



Silicate (Olivine stoichiometry) Amorphous Homogeneous 0.1, 1.5, 6.0 

Silicate (Pyroxene stoichiomtry) Amorphous Homogeneous 0.1, 1.5, 6.0 

Forsterite Crystalline Inhomogcncous 0.1, 1.5, 6.0 

Enstatite Crystalline Inhomogeneous 0.1, 1.5, 6.0 

Silica Amorphous Inhomogcncous 0.1, 1.5, 6.0 



Dorschner et al. (1995) 
Dorschner et al. (1995) 
Servoin & Piriou (1973) 
Jager et al. (1998) 
Spitzcr & Klcinman (1960) 



Note. — Dust species used to reproduce the 10 fim silicate feature. See Bouwman et al. (2006) for more details. 
All components were included in three different sizes. The shape refers to homogeneous spheres (with opacities 
derived from Mic theory) and to inhomogeneous spheres (using the distribution of hollow spheres from Min et al. 
2005). 



Table 4. Summary of Dust Properties from Model Fitting 



ID 






M /os 




M sll 




M t (M Moo „) 


T CO rit T ^ grain 




Soli 




Sfos 


CI 

^enst 




^amo 




Comments 



73-758 


0. 


O4-0-24 
y ^+0.04 


0. 


ru-o- 04 

u ^+0.23 


nnq- 004 
u.uuy +0 004 





01 ~ 001 

ui +0.01 





0. 


no-0.01 
yo +0.01 


0007~ a0001 

u - uuu ' +0.00002 


721 ~ 139 

' z± +552 


49i;f 




5. 


oo— 0.21 
zo +0.20 




— 


n no-0.03 
u - io +0.40 


1 


oc-0.95 
• zu +0.24 





5. 


91 —0.22 
zi +0.23 


62~ - 45 

u - uz, +0.38 






11-2146 


0. 


opr-0.29 
OJ +0.55 


0. 


uu +0.29 


02~ 001 

u - uz +0.01 





P19— 0.02 
• UZ +0.02 


01 _0 - 01 

u - ui +0.02 





ncr-0.02 
yJ +0.02 


n 9fi -0 - 22 

u - zu +0.49 


1286 + 324 


z ' d +186 




1 


07- 1.78 
'°'+2.81 


2. 


62" 2 ' 08 

uz +2.50 


n 4 f ?-°- 31 

u -^°+0.48 


1 


4 o-0.87 
•^°+0.06 


1 41-0.68 

i -^ i +0.09 


2. 


4 o-1.87 
•^°+1.90 


1 1 ^-o- 16 

± - id +0.15 






11-2037 


0. 


79— 0.36 
• ' z +0.18 


0. 


1O-017 
iy +0.34 


n n9-o oo 5 

u - uz +0.009 





pr-0.01 
• u,J +0.02 


09~ a01 

u - uz, +0.03 





Ol -0 04 
y± +0.02 


02~ 01 

u - uz, +0.06 


701 -93 
' oi +217 


ocr9-121 
odz +134 




5. 


90 — 1.44 
' ZO +0.41 


0. 


OO-0.27 
°°+3.63 


— 


1 


'^+0.06 


0-89+o:45 


4. 


9^-2.22 
zu +1.18 


1 OO" 016 

-■-•^+0.09 






11-2031 


0. 


70— 0.18 
• ' °+0.11 


0. 


9n-0.ro 

zu +0.15 


o^~ - 02 

u - UJ +0.02 





ni -0.01 
ul +0.01 


01 "° 01 
u - Ui +o.oi 


0. 


q4-0.02 
y ^+0.02 


02" 01 

u - uz +0.03 


723^ 


334;™ 




3 


• 08 +o:53 


0. 


1 1 -001 

■ ± - L +0.51 


1 07" - 18 

i - u '+0.10 





7C--0.64 
• ' J +0.63 


1 4n-°- 68 

i -^ u +0.10 


2 


p-i -1.02 
' d± +0.64 


1-09;^ 






14-183 


0. 


7q— 0.29 
• ' y +0.07 


0. 


nq-0.03 
UO +0.87 


u - uo +0.02 





1 9-O.O5 
,iz +0.06 


— 


0. 


09-O.O7 
' oz +0.05 


0-002;H° 1 


630; 29 


285 + 4i 




0. 


1 n-0.00 

±u +0.00 


0. 


1 n-0.00 

lu +0.00 


— 





1 1 —0.01 

• iJ -+o.ii 


— 


0. 


1 n-0.00 

• lu +o.oo 


1 1 _0 - 01 


SB1 




82-272 


0. 


no— 0.03 
uo +0.13 


0. 


0-1 —0.13 
o± +0.05 


06" - 01 
u - uu +o.oi 





07—0.02 
u '+0.02 


n no-0.02 
u - uo +0.02 





04-O.O4 
°^+0.03 


n nn-0.05 
u - uy +0.06 




i98;g 




1 


1 — 1.04 
• 10 +4.14 


0. 


^ u +0.45 


n 14-O.O4 


1 


01 —0.40 
• Oi +0.17 


n co-0.48 
U-O-J+0.45 


0. 


4fi-°-30 
^ u +0.45 


u - '°+0.12 


SB2 




12-2113 


0. 


_ n on 

ni _u -30 
• U1 +0.14 


0. 


_ , n 1 a 

07— U. 14 

°'+0.31 


n ni -o oi 

u - Ui +o.oi 





ni-0.01 
ui +0.03 


0.004;0;°? 3 


0. 


^ ^ n no 

no-0.03 

yo +o.oi 


u - uo +0.21 


628+ii 9 


291^6 




1 


Ofi-1.43 
•00+1.25 


1. 


14-0.97 


14-0- 04 

u - i ^+1.32 





n r-0.85 
yj +0.53 


1 oo -0 - 80 


1. 


fi7 -1.36 
' u '+1.07 


TO -0 39 

U - OZ +0.67 


SB1 




13-157 


0. 


QQ-0.00 

• uu +o.oo 


0. 


qo-0.02 
yo +0.01 


n nn4-o- 004 

u - uu ^+0.006 





01 -o.oi 

ui +0.02 


OO" 00 

u - uu +0.01 





no-0.02 
yo +0.01 


3fi-°- 04 

u - ou +0.06 


1 991 ~ 138 
izy± +153 


151+1 








0. 


n-0.00 
lu +0.00 


in- 0003 

u - iU +o.io 





i n-0.00 

• iU +o.oo 


n 9Q— 0.18 
U - ZO +0.99 


0. 


1 n-0.00 
lu +0.00 


u - id +0.80 






91-155 


0. 


1R-0.07 
lu +0.18 


0. 


81 -0- 18 
o± +0.07 


00" 00 
u - uu +o.oo 





no-O.Ol 

• uz +o.oi 


002" 002 
u.uuz +0 002 





07-0.01 

■ y ' +0.01 


n^-o- 04 

u - ud +2.75 


944- 80 , 


262;! 




0. 


O7-0.71 
y '+1.26 


5. 


r 4 -1.85 
d ^+0.42 




1 


19-0.56 
•^ z +0.07 


u -°'+0.57 


4. 


yi +0.29 


1 Q1-0.41 
i - zi +0.17 






54-1547 


0. 


74-0-10 
' ' ^+0.08 


0. 


14-0.10 

i4 +0.14 


OS" 02 

u - UJ +0.02 





nr-0.03 
■ UJ +0.04 


02" a01 

u - uz +0.02 


0. 


OO-0.05 
' oo +0.04 


0002" 001 

U.UUUZ_|_Q QQQ3 


091 -220 

yzi +350 


°^ u +110 




0. 


4O-0.38 
•^ y +1.91 


0. 


4 o-0.37 
48 +3.02 


°-28+ :34 


1 


19 -0.67 
• iz +0.33 


-■-•^'+0.03 


0. 


^4-0-42 
°^+2.01 


0.82~ ; 21 






13-1250 


0. 


c-o-0.48 
JO +0.33 


0. 


00-0.31 
°°+0.46 


02-° 01 
u - uz +o.oi 





n fi -0.02 
• uu +0.02 


006-°- 001 

U.UUO +0 006 


0. 


Ol -0.03 
yi +0.03 




790-144 
' zo +395 


o fi r-180 
oud +130 




2. 


co-2.47 
' 0,3 +1.55 


0. 


RO-0.60 
Dy +1.67 


21 -0 - 11 

u - zi +0.42 





OR" - 77 
■ y °+0.44 


u -° ' +1.13 


2. 


oo— 2.31 
°°+1.64 


n 77-0.56 
u - ' ' +0.34 






23-162 


0. 


no-0.08 
UO +0.34 


0. 


OQ-0.33 
oy +0.08 


o.oos;™ 03 


0.005:3:8°* 


02" - 01 


0. 


Q7-0.02 
y '+0.01 


1 1 0.09 


1010;^ 


1 9fF 38 
iyu +131 



Table 4 — Continued 



ID 






M fos 




M sil 




M t (M Moon ) 


T 

- 1 - cont 


T 

gram 






l~>py r 




C 




^amo 


^cry 


Comments 






i 03 -0 - 91 

i - u,3 +3.94 


J-- 1 ' +3.83 




U - 4O +0.91 


0.62J; 4 ? 


l.UO +30 4 


U - Od +0.52 






01-580 
DG-481 


u - oq: +0.08 

1 4" - 04 

u ' i ^+1.32 
3fi -011 

II- - 01 

u - iJ -+0.23 


49-°- 09 

u - 4Z +0.06 

o 1 1 -°- 01 

u - do +0.13 

o 4 c r°- 34 

U - 40 +3.13 


f) no" ' 005 
u - uz +o.oi 

2fi-°- 12 

u - zu +0.28 

n<r - 02 

J-- 1 '+0.25 


u - uz +0.01 
fl qq— 0.18 
U - ZO +0.84 

oi -° 01 

u - ui +0.03 

QS -0 - 84 
u - yo +0.47 


o.oo 2 ;H? 2 

i O0-1.49 

J^-^+o.ig 

oi -°- 01 

u - ui +0.02 

1 QQ-0.94 
i - oo +0.13 


Qfi" 02 
u - yu +0.01 

o n-° 03 

u - ±o +0.89 

n 04-0.05 

u - y ^+0.03 

0.23;?;}? 


oo pp - 001 

u.uuo +0 001 
n qq-o.18 

U.OO +0 45 

n<r - 01 

u - uo +0.01 

1 90-O.37 


1529;f 
13607™ 


233 +19 
1 77~ 9 

1 ' '+0 



Note. - Parameters and properties of the silicate from the fitting of the 10 /im feature (using the models 
described by Bouwman et al. 2006). Only the stars with good fits are included in the list. The parameters 
are the fraction of mass in the different species (olivine, pyroxene, forsterite, enstatite, silica, total amorphous 
grains, in units of fraction of the total mass), total mass in warm silicates (M t , in MMoon), continuum and grain 
temperature {T cont, T grain, in K), and mass-weighted average size of the grains of the different species (olivine, 
pyroxene, forsterite, enstatite, silica, total amorphous grains, and total crystalline grains in /im). All the quantities 
include their 1-a errors. 



Table 5. Correlations between stellar, disk and silicate characteristics 



Parameters 


Figure 


N 


r(low-mass) 


r(all) 


Possible Implications/Comments 


F8.6/F9.8 vs. F11.3/F9.8 


12 


20 


0.67 


0.67 


Similar to Bouwman et al. (2006) 


F8.6/F9.8 vs. F8.6/F11.3 




20 


0.37 


0.37 




F8.6/F9.8 vs. Spec. Type 




20 


0.40 


0.21 


Flux ratios correlation does not seem due to spectral type differences 


F8.6/F11.3 vs. Spec. Type 




20 


0.44 


0.25 




F11.3/F9.8 vs. Spec. Type 




20 


0.04 


0.02 




Fs.e/Fg.s vs. log(M) 




17 


0.40 


0.28 


Flux ratios are not significantly affected by the accretion rates 


Fs.e/Fii.3 vs. log(M) 




17 


0.35 


0.30 




F11.3/F9.8 vs. log(M) 




17 


0.25 


0.13 


55 


F8.6/F9.8 vs. F peafe 




20 


0.77 


0.52 


Weak correlation between the peak flux and silicate feature shape 


F8.6/F11.3 vs. F peak 




20 


0.52 


0.52 


55 


F11.3/F9.8 vs. F peak 


m 


20 


0.36 


0.36 


55 


Fp e<I fc vs. Spec. Type 


m 


20 


0.61 


0.22 


Higher peak in earlier stars (illumination? disk area?) but DG-481 is out of trend 


Fpeak VS M 




17 


0.13 


0.15 


The accretion rate does not affect the silicate strength. 


Fpeak vs log(M) 




17 


0.07 


0.17 


55 


Fpeak VS. H-K 




20 


0.04 


0.15 


No correlation peak flux - inner disk brightness | 


Fpeak VS H-[8] 




19 


0.22 


0.26 


The dependence of the silicate peak flux with the SED shape is very weak or noneCn 

55 

1 


Fpeak vs. H-[24] 


13 


20 


0.16 


0.14 


Fpeak VS H-[31] 




20 


0.35 


0.32 


55 


Fpeak vs [8]-[24] 




19 


0.24 


0.30 


55 


Fpeak vs [24]- [31] 


m 


19 


0.12 


0.15 




Fpeak vs. Age 


nu 


20 


0.46 


0.29 


Weak correlation peak flux - age for low-mass only 


Mass (total) vs. Spect. Type 


na 


13 


0.05 


0.13 


Total silicate mass does not depend on spectral type 


Mass (total) vs. Age 


ed 


13 


0.12 


0.10 


Total silicate mass does not depend significantly on age 


Log. Mass (total) vs. log( M) 


[161 


12 


0.78 


0.45 


More mass in silicates if M is larger (more turbulence?) 


Mass (amorphous) vs. H-K 




13 


0.62 


0.60 


Brighter inner disk have more mass in silicates 


Mass (total) vs. H-K 


M 


13 


0.62 


0.57 


55 


Mass (amorphous) vs. Stellar Mass 




13 


0.15 


0.22 


No correlation between mass in silicates and stellar mass 


Mass (total) vs. Stellar Mass 




13 


0.13 


0.16 




Mass (total) vs. M 




10 


0.71 


0.71 


More mass in silicates if M is larger (more turbulence?); (1) 


Mass (total) vs. log( M) 




12 


0.54 


0.43 




Mass (total) vs. Slope (H-[24]) 




12 


0.13 


0.18 


No correlation silicate/SED shape (slope) 


Size vs. Spec. Type 


M 


13 


0.39 


0.41 




Size (amorphous) vs. Age 


M 


13 


0.48 


0.52 


Large grains in atmosphere only at young ages (sedimentation?) 


Size (crystalline) vs. Age 


m 


13 


0.13 


0.12 


No correlation between the size of crystalline grains and the age 


Size vs. log( M) 


m 


12(11) 


0.39(15) 


0.11(0.35) 


(2) 



Table 5 — Continued 



Parameters 


Figure 


N 


r(low-mass) 


r(all) 


Possible Implications/Comments 


OlZC VS. 1V1 








12 


0.22 


0.15 




Size vs. F peak 






EE! 


13 


0.50 


0.38 


Smaller grains for larger peak fluxes (larger illumination or more flared disks?) 


Size (amorphous) 


vs. 


H-K 


1151 


1.3 


0.24 


0.11 


No correlation inner disk brightness/silicate size 


Size (amorphous) 


vs. 


H-[8] 


[15] 


1.3 


0.47 


0.48 


Smaller grains in more flared disks, sedimentation? 


Size (amorphous) 


vs. 


H-[24] 




1.3 


0.48 


0.49 




Size (amorphous) 


vs. 


[8]-[24] 


EES] 


1.3 


0.27 


0.12 




Size (amorphous) 


vs. 


[24]-[31] 


US] 


1.3 


0.48 


0.48 




Crystal. Fraction 


vs. 


Spec. Type 


[13 


1.3 


0.50 


0.23 


No clear dependence of crystallization on spectral type 


Crystal. Fraction 


vs. 


Age 


CE3 


1.3 


0.27 


0.19 


Crystallization does not change with age. Early formation? 


Crystal. Fraction 


vs. 


log( M) 


E] 


1.3 


0.18 


0.03 


Crystallinity fraction does not correlate with other properties 


Crystal. Fraction 


vs. 


H-K 


El 


13 


0.05 


0.06 




M vs. Age 








11(12) 


0.53(0.46) 


0.41(0.19) 


Accretion rate tends to decrease with time but there is large scatter (3) 


log( M) vs. Age 








11(12) 


0.46(0.18) 


0.19(0.20) 





Note. — Possible correlations and lacks of correlation between stellar, disk, and silicate properties investigated. If the correlation is displayed in any of the figures, 
the figure number is given. The significance of each checked relation is given by the correlation coefficient r, (including only low-mass stars, 'r(low-mass)', or low- 
and high-mass stars, 'r(all)'; note that the only high-mass star with measurable silicate is DG-481). Values of r = 1 would indicate a perfect correlation, and r = 
characterizes completely independent quantities. Note that all these relations are valid only within the limit of our data (for Class II objects with or without inner 
gaps, moderate to very low accretion rates, ages 1-12 Myr). The number of data points (N) includes only low-mass stars, add 1 if DG-481 is included. The accretion 
rates include upper limits (see Paper IV), and we assume M=10 — 12 MQ/yr for DG-481 (given the lack of Ho emission (Paper II). Note that 'Mass' correspond to the 
total mass in silicates in the warm disk atmosphere (given that only the warm dust gas in the upper disk layers produces the silicate feature). For all the stars, the mass 
fraction of crystalline silicates represents only a small fraction of the total mass; therefore we consider mostly the total mass. Similarly, we do not see any correlation 
between the sizes of crystalline grains and any other stellar/disk properties, so the parameter 'Size' refers only to mass-weigted amorphous grain size. (l)=Star 13-157 
not included; (2)= Star 73-758 not included; (3)= Values in parenthesis exclude G9 star 82-272, since G stars have larger accretion rates than K and M. 



Table 6. Excesses over the Photosphere of High- and Intermediate-Mass Stars 



ID 


kj^J^*.. . _L y 




lfWFin /Fin ^ 
lu 6^ r 10/ r lOp J 


lu 6^ r 24/ r 24p 7 


lnO"( Fon /Fin 1 

iu &V r 30/ r 30p ) 


iU &V r 70/ r 70pJ 




MVA-1312 


B4 


4 


0.0 


1.3 


2.1 





D 


KUN-314S 


A 


4 


1.0 


1.4 


1.4 





P: 


CCDM+5734A 


B3+B5 


4 


0.3 


1.0 


1.1 


— 


P 


MVA-468 


B7 


4 


0.0 


0.7 


1.0 




D 


MVA-426 


B7 


4 


2.0 


3.3 


3.9 


4.1 


P 


MVA-447 


FO 


4 


0.0 


1.0 


1.1 




D 


BD+572356 


A4 


4 


0.1 


0.4 


0.9 




D 


KUN-196 


B9 


4 


0.7 


1.4 


2.3 




D 


DG-481 


A7 


12 


0.7 


1.0 


1.2 




P 


DG-39 


AO 


12 


0.0 


0.9 


1.1 




D 


DG-682 


A2 


12 


0.0 


1.4 


1.4 




D 


DG-62 


F5 


12 












DG-912 


F5.5 


12 


0.0 


0.6 


1.3 




D 



Note. — Flux excess over the photosphere at different wavelengths (log(F;\/FAp, where F\ is the measured flux and F\ p is the 
photosphcric flux at the given wavelength) for the high- and intermediate-mass with disks. Note that the ages given are the average 
ages in the cluster. Errors arc 10%. No values are given for DG-62 due to the uncertainty of the measure and the possibility that 
the excess is not real. The disk type denotes debris- like disks ("D") or protoplanetary-likc disks ("P") based on the presence of 
silicate emission and in the magnitude of the excess. The colon (":") denotes uncertain values. 



